Single-Cell-Based High-Throughput Ig and TCR Repertoire Sequencing Analysis in Rhesus Macaques.

Recent advancements in microfluidics and high-throughput sequencing technologies have enabled recovery of paired H and L chains of Igs and VDJ and VJ chains of TCRs from thousands of single cells simultaneously in humans and mice. Despite rhesus macaques being one of the most well-studied model organisms for the human adaptive immune response, high-throughput single-cell immune repertoire sequencing assays are not yet available due to the complexity of these polyclonal receptors. We used custom primers that capture all known rhesus macaque Ig and TCR isotypes and chains that are fully compatible with a commercial solution for single-cell immune repertoire profiling. Using these rhesus-specific assays, we sequenced Ig and TCR repertoires in >60,000 cells from cryopreserved rhesus PBMCs, splenocytes, and FACS-sorted B and T cells. We were able to recover every Ig isotype and TCR chain, measure clonal expansion in proliferating T cells, and pair Ig and TCR repertoires with gene expression profiles of the same single cells. Our results establish the ability to perform high-throughput immune repertoire analysis in rhesus macaques at the single-cell level.

Single-Cell Transcriptomic Profiling of MAIT Cells in Patients With COVID-19.

Background: Mucosal-associated invariant T (MAIT) cells are considered to participate of the host immune response against acute severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection; however, single-cell transcriptomic profiling of MAIT cells in patients with COVID-19 remains unexplored. Methods: We performed single-cell RNA sequencing analyses on peripheral MAIT cells from 13 patients with COVID-19 and 5 healthy donors. The transcriptional profiles of MAIT cells, together with assembled T-cell receptor sequences, were analyzed. Flow cytometry analysis was also performed to investigate the properties of MAIT cells. Results: We identified that differentially expressed genes (DEGs) of MAIT cells were involved in myeloid leukocyte activation and lymphocyte activation in patients with COVID-19. In addition, in MAIT cells from severe cases, more DEGs were enriched in adaptive cellular and humoral immune responses compared with those in moderate cases. Further analysis indicated that the increase of cell cytotoxicity (killing), chemotaxis, and apoptosis levels in MAIT cells were consistent with disease severity and displayed the highest levels in patients with severe disease. Interestingly, flow cytometry analysis showed that the frequencies of pyroptotic MAIT cells, but not the frequencies of apoptotic MAIT cells, were increased significantly in patients with COVID-19, suggesting pyroptosis is one of leading causes of MAIT cell deaths during SARS-CoV-2 infection. Importantly, there were more clonal expansions of MAIT cells in severe cases than in moderate cases. Conclusions: The results of the present study suggest that MAIT cells are likely to be involved in the host immune response against SARS-CoV-2 infection. Simultaneously, the transcriptomic data from MAIT cells provides a deeper understanding of the immune pathogenesis of the disease.

Contribution of T Cell Receptor Alpha and Beta CDR3, MHC Typing, V and J Genes to Peptide Binding Prediction.

Introduction: Predicting the binding specificity of T Cell Receptors (TCR) to MHC-peptide complexes (pMHCs) is essential for the development of repertoire-based biomarkers. This affinity may be affected by different components of the TCR, the peptide, and the MHC allele. Historically, the main element used in TCR-peptide binding prediction was the Complementarity Determining Region 3 (CDR3) of the beta chain. However, recently the contribution of other components, such as the alpha chain and the other V gene CDRs has been suggested. We use a highly accurate novel deep learning-based TCR-peptide binding predictor to assess the contribution of each component to the binding. Methods: We have previously developed ERGO-I (pEptide tcR matchinG predictiOn), a sequence-based T-cell receptor (TCR)-peptide binding predictor that employs natural language processing (NLP) -based methods. We improved it to create ERGO-II by adding the CDR3 alpha segment, the MHC typing, V and J genes, and T cell type (CD4+ or CD8+) as to the predictor. We then estimate the contribution of each component to the prediction. Results and Discussion: ERGO-II provides for the first time high accuracy prediction of TCR-peptide for previously unseen peptides. For most tested peptides and all measures of binding prediction accuracy, the main contribution was from the beta chain CDR3 sequence, followed by the beta chain V and J and the alpha chain, in that order. The MHC allele was the least contributing component. ERGO-II is accessible as a webserver at http://tcr2.cs.biu.ac.il/ and as a standalone code at https://github.com/IdoSpringer/ERGO-II.

Global characterization of B cell receptor repertoire in COVID-19 patients by single-cell V(D)J sequencing.

The world is facing a pandemic of Corona Virus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Adaptive immune responses are essential for SARS-CoV-2 virus clearance. Although a large body of studies have been conducted to investigate the immune mechanism in COVID-19 patients, we still lack a comprehensive understanding of the BCR repertoire in patients. In this study, we used the single-cell V(D)J sequencing to characterize the BCR repertoire across convalescent COVID-19 patients. We observed that the BCR diversity was significantly reduced in disease compared with healthy controls. And BCRs tend to skew toward different V gene segments in COVID-19 and healthy controls. The CDR3 sequences of heavy chain in clonal BCRs in patients were more convergent than that in healthy controls. In addition, we discovered increased IgG and IgA isotypes in the disease, including IgG1, IgG3 and IgA1. In all clonal BCRs, IgG isotypes had the most frequent class switch recombination events and the highest somatic hypermutation rate, especially IgG3. Moreover, we found that an IgG3 cluster from different clonal groups had the same IGHV, IGHJ and CDR3 sequences (IGHV4-4-CARLANTNQFYDSSSYLNAMDVW-IGHJ6). Overall, our study provides a comprehensive characterization of the BCR repertoire in COVID-19 patients, which contributes to the understanding of the mechanism for the immune response to SARS-CoV-2 infection.

Chromosomal translocation t(11;14) and p53 deletion induced by the CRISPR/Cas9 system in normal B cell-derived iPS cells.

Multiple myeloma (MM) cells are derived from mature B cells based on immunoglobulin heavy chain (IgH) gene analysis. The onset of MM is often caused by a reciprocal chromosomal translocation (cTr) between chr 14 with IgH and chr 11 with CCND1. We propose that mature B cells gain potential to transform by reprograming, and then chromosomal aberrations cause the development of abnormal B cells as a myeloma-initiating cell during B cell redifferentiation. To study myeloma-initiating cells, we have already established normal B cell-derived induced pluripotent stem cells (BiPSCs). Here we established two BiPSCs with reciprocal cTr t(11;14) using the CRISPR/Cas9 system; the cleavage site were located in the IgH Eµ region of either the VDJ rearranged allele or non-rearranged allele of IgH and the 5'-upsteam region of the CCND1 (two types of BiPSC13 with t(11;14) and MIB2-6 with t(11;14)). Furthermore, p53 was deleted using the CRISPR/Cas9 system in BiPSC13 with t(11;14). These BiPSCs differentiated into hematopoietic progenitor cells (HPCs). However, unlike cord blood, those HPCs did not differentiated into B lymphocytes by co-culture with BM stromal cell. Therefore, further ingenuity is required to differentiate those BiPSCs-derived HPCs into B lymphocytes.

The T Cell Repertoires from Nickel Sensitized Joint Implant Failure Patients.

Nickel (Ni2+) is one of the most common allergens, affecting around 10-15% of the general population. As the demand for orthopedic implant surgery rises, the number of surgical revisions due to joint implant failure also increases. There is evidence that some patients develop joint failure due to an immune response to a component of the implant, and we have found that Ni2+ is an especially important cause. Hence, understanding the mechanisms by which Ni2+ allergy induces joint implant failure becomes a critical research question. The structural basis of Ni2+ activation of pathogenic T cells is still not clear. The purpose of this study was to characterize Ni2+-reactive T cell repertoires derived from the peripheral blood of joint failure patients due to Ni2+ sensitization using single-cell sequencing techniques. We stimulated the proliferation of Ni2+ -reactive T cells from two implant failure patients in vitro, and sorted them for single-cell VDJ sequencing (10× genomics). We identified 2650 productive V-J spanning pairs. Both TCR α chains and ß chains were enriched. TRBV18 usage is the highest in the P7 CD4+ population (18.1%), and TRBV5-1 usage is the highest in the P7 CD8+ population (12.1%). TRBV19 and TRBV20-1 segments are present in a high percentage of both P7 and P9 sequenced T cells. Remarkably, the alpha and beta chain combination of TRAV41-TRBV18 accounts for 13.5% of the CD4+ population of P7 patient. Compared to current Ni specific T cell repertoire studies of contact dermatitis, the Vα and Vß usages of these joint implant failure patients were different. This could be due to the different availability of self-peptides in these two different tissues. However, TRBV19 (Vß17) was among frequently used TCR ß chains, which are common in previous reports. This implies that some pathogenic T cells could be similar in Ni2+ hypersensitivities in skin and joints. The alignment of the TCR CDR3ß sequences showed a conserved glutamic acid (Glu) that could potentially interact with Ni2+. The study of these Ni2+ specific TCRs may shed light on the molecular mechanism of T cell activation by low molecular weight chemical haptens.

T Cell Receptor Beta-Chain Profiling of Tumor Tissue, Peripheral Blood and Regional Lymph Nodes From Patients With Papillary Thyroid Carcinoma.

Objective: To study the characteristics of the T cell receptor (TCR) repertoire in cancer tissue, peripheral blood and regional lymph nodes (LNs) from patients with papillary thyroid carcinoma (PTC). Methods: PTC tissue, peripheral blood mononuclear cells (PBMCs) and regional LNs of six patients with papillary thyroid carcinoma were harvested. T cell receptor beta-chain (TCRß) profiling was performed though high-throughput sequencing (HTS), and IMonitor, MiXCR and VDJtools were used to analyze the characteristics of the TCR repertoire. Results: The results of IMonitor and those of MiXCR and VDJtools were very similar. The unique CDR3 of TCRß from LNs was higher than that of PBMCs, and the CDR3 of TCRß from LNs was higher than that of PTC tissue. Shannon's diversity index, D50, inverse Simpson index_mean and normalized Shannon's diversity index_mean of CDR3 from LNs were higher than those of PTCs and PBMCs. The HEC (high expansion clones) rate of CDR3 sequences at the amino acid level in PTC tissue was higher than that of PBMCs, which was higher than that of LNs. The V-J HEC rate of CDR3 was highest in PTC tissue, followed by PBMCs and LNs. Conclusion: TCR CDR3 profiling showed differences among and within the PBMCs, PTC tissues and regional LNs of PTC, including unique CDR3, CDR3 HEC at the amino acid level, CDR3 V-J HEC at the amino acid level, Shannon's diversity index and D50. The TCRß repertoire of PTC tissue, peripheral blood and regional LNs of PTC provide a reference for further study of immunity mechanisms against PTC.

The thymoproteasome hardwires the TCR repertoire of CD8+ T cells in the cortex independent of negative selection.

The thymoproteasome expressed specifically in thymic cortical epithelium optimizes the generation of CD8+ T cells; however, how the thymoproteasome contributes to CD8+ T cell development is unclear. Here, we show that the thymoproteasome shapes the TCR repertoire directly in cortical thymocytes before migration to the thymic medulla. We further show that the thymoproteasome optimizes CD8+ T cell production independent of the thymic medulla; independent of additional antigen-presenting cells, including medullary thymic epithelial cells and dendritic cells; and independent of apoptosis-mediated negative selection. These results indicate that the thymoproteasome hardwires the TCR repertoire of CD8+ T cells with cortical positive selection independent of negative selection in the thymus.

Alignment free identification of clones in B cell receptor repertoires.

Following antigenic challenge, activated B cells rapidly expand and undergo somatic hypermutation, yielding groups of clonally related B cells with diversified immunoglobulin receptors. Inference of clonal relationships based on the receptor sequence is an essential step in many adaptive immune receptor repertoire sequencing studies. These relationships are typically identified by a multi-step process that involves: (i) grouping sequences based on shared V and J gene assignments, and junction lengths and (ii) clustering these sequences using a junction-based distance. However, this approach is sensitive to the initial gene assignments, which are error-prone, and fails to identify clonal relatives whose junction length has changed through accumulation of indels. Through defining a translation-invariant feature space in which we cluster the sequences, we develop an alignment free clonal identification method that does not require gene assignments and is not restricted to a fixed junction length. This alignment free approach has higher sensitivity compared to a typical junction-based distance method without loss of specificity and PPV. While the alignment free procedure identifies clones that are broadly consistent with the junction-based distance method, it also identifies clones with characteristics (multiple V or J gene assignments or junction lengths) that are not detectable with the junction-based distance method.

Potent Neutralizing Antibodies against SARS-CoV-2 Identified by High-Throughput Single-Cell Sequencing of Convalescent Patients' B Cells.

The COVID-19 pandemic urgently needs therapeutic and prophylactic interventions. Here, we report the rapid identification of SARS-CoV-2-neutralizing antibodies by high-throughput single-cell RNA and VDJ sequencing of antigen-enriched B cells from 60 convalescent patients. From 8,558 antigen-binding IgG1+ clonotypes, 14 potent neutralizing antibodies were identified, with the most potent one, BD-368-2, exhibiting an IC50 of 1.2 and 15 ng/mL against pseudotyped and authentic SARS-CoV-2, respectively. BD-368-2 also displayed strong therapeutic and prophylactic efficacy in SARS-CoV-2-infected hACE2-transgenic mice. Additionally, the 3.8 Å cryo-EM structure of a neutralizing antibody in complex with the spike-ectodomain trimer revealed the antibody's epitope overlaps with the ACE2 binding site. Moreover, we demonstrated that SARS-CoV-2-neutralizing antibodies could be directly selected based on similarities of their predicted CDR3H structures to those of SARS-CoV-neutralizing antibodies. Altogether, we showed that human neutralizing antibodies could be efficiently discovered by high-throughput single B cell sequencing in response to pandemic infectious diseases.

Comparisons of the antibody repertoires of a humanized rodent and humans by high throughput sequencing.

The humanization of animal model immune systems by genetic engineering has shown great promise for antibody discovery, tolerance studies and for the evaluation of vaccines. Assessment of the baseline antibody repertoires of unimmunized model animals will be useful as a benchmark for future immunization experiments. We characterized the heavy chain and kappa light chain antibody repertoires of a model animal, the OmniRat, by high throughput antibody sequencing and made use of two novel datasets for comparison to human repertoires. Intra-animal and inter-animal repertoire comparisons reveal a high level of conservation in antibody diversity between the lymph node and spleen and between members of the species. Multiple differences were found in both the heavy and kappa chain repertoires between OmniRats and humans including gene segment usage, CDR3 length distributions, class switch recombination, somatic hypermutation levels and in features of V(D)J recombination. The Inference and Generation of Repertoires (IGoR) software tool was used to model recombination in VH regions which allowed for the quantification of some of these differences. Diversity estimates of the OmniRat heavy chain repertoires almost reached that of humans, around two orders of magnitude less. Despite variation between the species repertoires, a high frequency of OmniRat clonotypes were also found in the human repertoire. These data give insights into the development and selection of humanized animal antibodies and provide actionable information for use in vaccine studies.

Concurrent immunoglobulin G-lambda type multiple myeloma and mixed cellularity classical Hodgkin lymphoma: A case report.

A 66-year-old man with a swollen right inguinal lymph node (LN) had pain on the lower side of the back. Computed tomography revealed bone disease in the back and swollen right inguinal LNs. Laboratory studies showed anemia and serum immunoglobulin G-lambda (IgG-λ) type monoclonal protein. The bone marrow contained 39.6% plasma cells. He was diagnosed with IgG-λ type multiple myeloma (MM). However, the pathological findings of the right inguinal LN were mixed cellular classical Hodgkin lymphoma (HL). The administration of melphalan, prednisone, and bortezomib (MPB) was started for MM; however, swelling in the right inguinal LN increased. After three cycles of MPB, the administration of doxorubicin, bleomycin, vinblastine, and dacarbazine (ABVD) was started for HL. However, HL was refractory to ABVD. Pancytopenia subsequently progressed and rapid swelling occurred in his LNs. He died 7 months after diagnosis. Multiple myeloma was diagnosed, based on the typical symptoms, although the pathological findings of the LN indicated a diagnosis of HL. We analyzed the molecular relationship between MM and HL cells using a direct sequencing method. The sequencing results demonstrated that the variable-diversity-joining (VDJ) region of the IgH gene was identified with 94.4% of IGLV3-32*01 in the bone marrow sample at diagnosis. Furthermore, clonotypic IgH sequence was identified in CD30-positive cells from the LN. These results suggested that the clonal HL cells were derived from the same source as the clonal MM cells and demonstrated that MM and HL in this patient may have originated from the same B cell progenitor.

Characterization of novel, recurrent genomic rearrangements as sensitive MRD targets in childhood B-cell precursor ALL.

B-cell precursor (BCP) ALL carry a variety of classical V(D)J rearrangements as well as genomic fusions and translocations. Here, we assessed the value of genomic capture high-throughput sequencing (gc-HTS) in BCP ALL (n = 183) for the identification and implementation of targets for minimal residual disease (MRD) testing. For TRδ, a total of 300 clonal rearrangements were detected in 158 of 183 samples (86%). Beside clonal Vδ2-Dδ3, Dδ2-Dδ3, and Vδ2-Jα we identified a novel group of recurrent Dδ-Jα rearrangements, comprising Dδ2 or Dδ3 segments fused predominantly to Jα29. For IGH-JH, 329 clonal rearrangements were identified in 172 of 183 samples (94%) including novel types of V(D)J joining. Oligoclonality was found in ~1/3 (n = 57/183) of ALL samples. Genomic breakpoints were identified in 71 BCP-ALL. A distinct MRD high-risk subgroup of IGH-V(D)J-germline ALL revealed frequent deletions of IKZF1 (n = 7/11) and the presence of genomic fusions (n = 10/11). Quantitative measurement using genomic fusion breakpoints achieved equivalent results compared to conventional V(D)J-based MRD testing and could be advantageous upon persistence of a leukemic clone. Taken together, selective gc-HTS expands the spectrum of suitable MRD targets and allows for the identification of genomic fusions relevant to risk and treatment stratification in childhood ALL.

Stability and uniqueness of clonal immunoglobulin CDR3 sequences for MRD tracking in multiple myeloma.

Minimal residual disease (MRD) tracking, by next generation sequencing of immunoglobulin sequences, is moving towards clinical implementation in multiple myeloma. However, there is only sparse information available to address whether clonal sequences remain stable for tracking over time, and to what extent light chain sequences are sufficiently unique for tracking. Here, we analyzed immunoglobulin repertoires from 905 plasma cell myeloma and healthy control samples, focusing on the third complementarity determining region (CDR3). Clonal heavy and/or light chain expression was identified in all patients at baseline, with one or more subclones related to the main clone in 3.2%. In 45 patients with 101 sequential samples, the dominant clonal CDR3 sequences remained identical over time, despite differential clonal evolution by whole exome sequencing in 49% of patients. The low frequency of subclonal CDR3 variants, and absence of evolution over time in active multiple myeloma, indicates that tumor cells at this stage are not under selective pressure to undergo antibody affinity maturation. Next, we establish somatic hypermutation and non-templated insertions as the most important determinants of light chain clonal uniqueness, identifying a potentially trackable sequence in the majority of patients. Taken together, we show that dominant clonal sequences identified at baseline are reliable biomarkers for long-term tracking of the malignant clone, including both IGH and the majority of light chain clones.

Detection of a Low Level and Heterogeneous B Cell Immune Response in Peripheral Blood of Acute Borreliosis Patients With High Throughput Sequencing.

The molecular diagnosis of acute Borreliosis is complicated and better strategies to improve the diagnostic processes are warranted. High Throughput Sequencing (HTS) of human B cell repertoires after e.g., Dengue virus infection or influenza vaccination revealed antigen-associated "CDR3 signatures" which may have the potential to support diagnosis in infectious diseases. The human B cell immune response to Borrelia burgdorferi sensu lato-the causative agent of Borreliosis-has mainly been studied at the antibody level, while less attention has been given to the cellular part of the humoral immune response. There are indications that Borrelia actively influence the B cell immune response and that it is therefore not directly comparable to responses induced by other infections. The main goal of this study was to identify B cell features that could be used to support diagnosis of Borreliosis. Therefore, we characterized the B cell immune response in these patients by combining multicolor flow cytometry, single Borrelia-reactive B cell receptor (BCR) sequencing, and B cell repertoire deep sequencing. Our phenotyping experiments showed, that there is no significant difference between B cell subpopulations of acute Borreliosis patients and controls. BCR sequences from individual epitope-reactive B cells had little in common between each other. HTS showed, however, a higher complementarity determining region 3 (CDR3) amino acid (aa) sequence overlap between samples from different timepoints in patients as compared to controls. This indicates, that HTS is sensitive enough to detect ongoing B cell immune responses in these patients. Although each individual's repertoire was dominated by rather unique clones, clustering of bulk BCR repertoire sequences revealed a higher overlap of IgG BCR repertoire sequences between acute patients than controls. Even if we have identified a few Borrelia-associated CDR3aa sequences, they seem to be rather unique for each patient and therefore not suitable as biomarkers.

Inferred Allelic Variants of Immunoglobulin Receptor Genes: A System for Their Evaluation, Documentation, and Naming.

Immunoglobulins or antibodies are the main effector molecules of the B-cell lineage and are encoded by hundreds of variable (V), diversity (D), and joining (J) germline genes, which recombine to generate enormous IG diversity. Recently, high-throughput adaptive immune receptor repertoire sequencing (AIRR-seq) of recombined V-(D)-J genes has offered unprecedented insights into the dynamics of IG repertoires in health and disease. Faithful biological interpretation of AIRR-seq studies depends upon the annotation of raw AIRR-seq data, using reference germline gene databases to identify the germline genes within each rearrangement. Existing reference databases are incomplete, as shown by recent AIRR-seq studies that have inferred the existence of many previously unreported polymorphisms. Completing the documentation of genetic variation in germline gene databases is therefore of crucial importance. Lymphocyte receptor genes and alleles are currently assigned by the Immunoglobulins, T cell Receptors and Major Histocompatibility Nomenclature Subcommittee of the International Union of Immunological Societies (IUIS) and managed in IMGT®, the international ImMunoGeneTics information system® (IMGT). In 2017, the IMGT Group reached agreement with a group of AIRR-seq researchers on the principles of a streamlined process for identifying and naming inferred allelic sequences, for their incorporation into IMGT®. These researchers represented the AIRR Community, a network of over 300 researchers whose objective is to promote all aspects of immunoglobulin and T-cell receptor repertoire studies, including the standardization of experimental and computational aspects of AIRR-seq data generation and analysis. The Inferred Allele Review Committee (IARC) was established by the AIRR Community to devise policies, criteria, and procedures to perform this function. Formalized evaluations of novel inferred sequences have now begun and submissions are invited via a new dedicated portal (https://ogrdb.airr-community.org). Here, we summarize recommendations developed by the IARC-focusing, to begin with, on human IGHV genes-with the goal of facilitating the acceptance of inferred allelic variants of germline IGHV genes. We believe that this initiative will improve the quality of AIRR-seq studies by facilitating the description of human IG germline gene variation, and that in time, it will expand to the documentation of TR and IG genes in many vertebrate species.

Mutating for Good: DNA Damage Responses During Somatic Hypermutation.

Somatic hypermutation (SHM) of immunoglobulin (Ig) genes plays a key role in antibody mediated immunity. SHM in B cells provides the molecular basis for affinity maturation of antibodies. In this way SHM is key in optimizing antibody dependent immune responses. SHM is initiated by targeting the Activation-Induced Cytidine Deaminase (AID) to rearranged V(D)J and switch regions of Ig genes. The mutation rate of this programmed mutagenesis is ~10-3 base pairs per generation, a million-fold higher than the non-AID targeted genome of B cells. AID is a processive enzyme that binds single-stranded DNA and deaminates cytosines in DNA. Cytosine deamination generates highly mutagenic deoxy-uracil (U) in the DNA of both strands of the Ig loci. Mutagenic processing of the U by the DNA damage response generates the entire spectrum of base substitutions characterizing SHM at and around the initial U lesion. Starting from the U as a primary lesion, currently five mutagenic DNA damage response pathways have been identified in generating a well-defined SHM spectrum of C/G transitions, C/G transversions, and A/T mutations around this initial lesion. These pathways include (1) replication opposite template U generates transitions at C/G, (2) UNG2-dependent translesion synthesis (TLS) generates transversions at C/G, (3) a hybrid pathway comprising non-canonical mismatch repair (ncMMR) and UNG2-dependent TLS generates transversions at C/G, (4) ncMMR generates mutations at A/T, and (5) UNG2- and PCNA Ubiquitination (PCNA-Ub)-dependent mutations at A/T. Furthermore, specific strand-biases of SHM spectra arise as a consequence of a biased AID targeting, ncMMR, and anti-mutagenic repriming. Here, we review mammalian SHM with special focus on the mutagenic DNA damage response pathways involved in processing AID induced Us, the origin of characteristic strand biases, and relevance of the cell cycle.

Non-coding Class Switch Recombination-Related Transcription in Human Normal and Pathological Immune Responses.

Antibody class switch recombination (CSR) to IgG, IgA, or IgE is a hallmark of adaptive immunity, allowing antibody function diversification beyond IgM. CSR involves a deletion of the IgM/IgD constant region genes placing a new acceptor Constant gene, downstream of the VDJH exon. CSR depends on non-coding (CSRnc) transcription of donor Iµ and acceptor IH exons, located 5' upstream of each CH coding gene. Although, our knowledge of the role of CSRnc transcription has advanced greatly, its extension and importance in healthy and diseased humans is scarce. We analyzed CSRnc transcription in 70,603 publicly available RNA-seq samples, including GTEx, TCGA, and the Sequence Read Archive using recount2, an online resource consisting of normalized RNA-seq gene and exon counts, as well as, coverage BigWig files that can be programmatically accessed through R. CSRnc transcription was validated with a qRT-PCR assay for Iµ, Iγ3, and Iγ1 in humans in response to vaccination. We mapped IH transcription for the human IGH locus, including the less understood IGHD gene. CSRnc transcription was restricted to B cells and is widely distributed in normal adult tissues, but predominant in blood, spleen, MALT-containing tissues, visceral adipose tissue and some so-called "immune privileged" tissues. However, significant Iγ4 expression was found even in non-lymphoid fetal tissues. CSRnc expression in cancer tissues mimicked the expression of their normal counterparts, with notable pattern changes in some common cancer subsets. CSRnc transcription in tumors appears to result from tumor infiltration by B cells, since CSRnc transcription was not detected in corresponding tumor-derived immortal cell lines. Additionally, significantly increased Iδ transcription in ileal mucosa in Crohn's disease with ulceration was found. In conclusion, CSRnc transcription occurs in multiple anatomical locations beyond classical secondary lymphoid organs, representing a potentially useful marker of effector B cell responses in normal and pathological immune responses. The pattern of IH exon expression may reveal clues of the local immune response (i.e., cytokine milieu) in health and disease. This is a great example of how the public recount2 data can be used to further our understanding of transcription, including regions outside the known transcriptome.